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ABSTRACT 

o : 

^ ' We estimate the strength of the bandpass-integrated thermal emission from 

the extrasolar planet HD 209458b at 3.6, 4.5, 5.8, and 8.0 /mi using the In- 
frared Array Camera (IRAC) on the Spitzer Space Telescope. We observe a 
single secondary eclipse simultaneously in all four bandpasses and find relative 
eclipse depths of 0.00094 ± 0.00009, 0.00213 ± 0.00015, 0.00301 ± 0.00043, and 
0.00240 ± 0.00026, respectively. These eclipse depths reveal that the shape of 
the inferred emission spectrum for the planet differs significantly from the pre- 
dictions of standard atmosphere models; instead the most plausible explanation 
would require the presence of an inversion layer high in the atmosphere lead- 
ing to significant water emission in the 4.5 and 5.8 /xm bandpasses. This is the 
first clear indication of such a temperature inversion in the atmosphere of a hot 
Jupiter, as previous observations of other planets appeared to be in reasonably 
good agreement with the predictions of models without such an inversion layer. 
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1. Introduction 

"Hot Jupiters" are a class of gas giant planets orbiting extremely close to their par- 



ent stars. Twenty one of these planets 



rave been observed to transit t heir parent stars 



(McCullough et a 




2or 


6 


: O'Donovan et al. 


2006. 


2007 


Bakos et al. 


2007al 


D:ICollier Cameron et al. 


2007; 


Burke et al. 


2007 


: 


Charbonneau et al.ll2007a 


Mandushev et al. 


2007 


Torres et al.ll2007 





allowing us to estimate not only their masses but also radii, temperatures, and other prop- 
erties. Of these twenty one systems, fifteen are bright enough for the type of Spitzer ob- 
servations described in this paper, including: HAT-P-1,2, and 3, HD 189733, HD 209458, 
HD 149026, TrES-1, 2, 3, and 4, WASP-1 and 2, and XO-1, 2, and 3. Of the fainter 
systems, including the five OGLE planets and CoRoT-Exo-1, most if not all (the position 
of CoRoT-Exo-1 has not yet been announced) lie in crowded fields where it is possible to 
achieve high quality infrared relative photometry using la rge-aperture ground-base d tele- 
scopes and nearby comparison stars, as demonstrated by ISnellen fc Covind (120071 ). The 
radii of most of these planets are well-described by basic models of an irradiated hot Jupiter 



( Bodcnheimer et al 



2003 



2001 



Laughlin et al. 



2005a 



Showman fc Guillotll2002; Bodenheimer et al. 2003: Baraffe et al. 



Burrows et al.ll2007al ). but there are currently four planets (HD 



209458b, WASP -lb, TrES-2, and 



than predicted (lO'Donovan et al. 



iES- 4 ) with radii that appear to be significantly larger 



2006 



Knutson et al. 2007a; Charbonneau et al. 2007b 



Shporer et al.l 120071 ; iMandushev et al.l 120071 ). TrES-4 is the most extr eme example, with 



an av erage density of on 



y 0.222 ±0.045 g cm" 3 (IMandushev et alJl2007h . The planet HAT- 



P-lb (IBakos et al.ll2007al) was i nitially thought to share this property as well, but subsequent 



observations (I Winn et al.ll2007l ) revealed that its radius was smaller than initially estimated. 



Given that there is a clear distinction between hot Jupiters with inflated radii and 
those with normal (i.e. consistent with models) radii, a comparative study of the infrared 
emission spectra of the planets in these two classes might reveal important differences. There 
are currently published estimates of the infrared emission from five planets, although this 



been observed at 8.0 /im, 7 .5 



Charbonneau 


et al. 


2005 


). and HD 189733b has 


and 16 urn ( 


Deming et al. 


2006 




Grillmair et al. 



20071 ; iKnutson et al.ll2007bl ). HD 209458b, the only planet with an inflated radius in this 



samp le, has been observed at 7.5 — 13.2 /im and 24 /xm (IDeming et al 



20071 ). Lastly, the core-dominated hot Jupiter HD 149026b ([Harrington et al.l 120071 ) and the 



2005 



Richardson et al. 
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hot Neptune G J 436b (jDeming et al.l l2007bl ; iDemory et al.l 120071 ) have both been observed 
at 8.0 /xm. This last planet is of particular interest, both because of its small size and 
significantly lower 8.0 /im brightness temperature of ~ 710 K. It is likely that this planet has 
an atmospheric composition significantly different than its hotter and more massive cousins, 
which should be revealed by future Spitzer observations at additional wavelengths. 

Based on these initial results, it appears that the dayside emission spectra for the 



tent with the predictions of standard cloud-free atmosphere models ( 


Sudarskv et al. 2003; 


Seaeer et al. 2005: Barman et al. 2005; Fortnev et al. 2005; Burrows et al. 


2006 


). The single 


exception is the 7.5 — 14.7 um IRS spectrum of HD 189733b measured by 


Grillmair et al. 



(120071 ). which did not appe ar to have the predict ed water absorption feature at the short- 
est wavelengths. However, iFortney et al.l (12006a ) suggested that this planet might have a 
close-t o-isothermal temperature profile on the dayside, which would wash out any absorption 
bands. IFortney fc Marleyl (120071 ) note that the broadband 8 /im eclipse depth for this planet 
appears to be deeper than the comparable eclipse depth from the IRS spectrum at these wave- 
lengths, which suggests imperfect instrumental calibration of the IRS specra. Intriguingly, 
the 7.5 — 13.2 /im IRS spectrum of the inflated planet HD 209458b also shows no evidence of 
a water absorption feature at the shortest wavelengths; in addition, this spectrum appears 
to contain two emiss ion features, one of wh ich has been tentatively identified as emission 
from silicate clouds (IRichardson et al.l 120071 ). The presence of high-altitude clouds in the 
atmosphere of HD 209458b would also explain the weaker-than-pr edicted sodium absorption 
in the planet's transmission spectrum (jCharbonneau et al.l 120021 ). but there is currently no 
definitive evidence for the presence of such clouds. 

In this paper we present the first estimates of the infrared emission spectrum of a planet 
from the "inflated" category of hot Jupiters at wavelengths shorter than 7.5 /im. Atmosphere 
models for these planets predict that many of the strongest features from CO, H 2 0, and CH 4 
will be located at wavelengths shorter than 8.0 /im, making this a particularly interesting 
wavelength region to study. The higher precision achievable in the 3 — 8 /mi wavelength 
range and complementary mult i- wavelength information allows us to test the predictions of 
atmosphere models for this planet, including the pressure-temperature profile of the dayside 
atmosphere and the relative strength of potential emission and/or absorption features in the 
planet's spectrum from CO, H 2 0, and CH 4 . 
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2. Observations 

We observed HD 209458b over a period of 8.1 hours on UT 200 5 Nov. 27, 



a single secondary eclipse, us ing the Infrared Ar ray Camera (IRAC) (IFazio et al.ll2004l ) on 



spann ing 



the Spitzer Space Telescope (IWerner et al.l I2004J ) . We observed in subarray mode with an 



exposure time of 0.1 s and cycled between the four IRAC channels in order to obtain estimates 
of the depth of the eclipse at 3.6, 4.5, 5.8, and 8.0 /im simultaneously. We obtained a total 
of 35,840 32 x 32 pixel images in each channel]. Images are taken in sets of 64, and four sets 
of images (4 x 64 images total) are obtained in each channel before repointing the telescope 
to position the star correctly on the subarray for the next channel. The position of the star 
on the subarray is still varying through the first set of 64 images after each repointing, and 
we chose to discard this initial set of images, leaving a total of 26,880 usable images in each 
of the four channels. The total size of this pointing drift is 0.3 pixels in the first set of 64 
images, and 0.1 pixels or less in the following three sets of images. 

Because the two shortest wavelength IRAC channels (3.6 and 4.5 jttm) use InSb detectors 
and the two longer wavelength channels (5.8 and 8.0 /mi) use Si:As detectors, there are 
fundamental differences between the properties of the data taken with these two types of 
detectors. We describe our analysis for each type of detector separately below. 



2.1. 3.6 and 4.5 fira Observations (InSb Detector) 

Because HD 209458 is a bright star (K = 6.31) and the background at these shorter 
wavelengths is minimal, we calculate the flux from the star in each image using aperture 
photometry with a radius of five pixels. We determine the position of the star in each image 
as the position-weighted sum of the flux in a 7 x 7 pixel box centered on the approximate 
position of the star. We estimate the background in each image by selecting a subset of 
pixels from the corners of the image where the point spread function of the star is faintest, 
making a histogram of the flux values in these pixels, and fitting a Gaussian function to the 
center of this distribution. We calculate the JD value for each image as the time at mid- 
exposure, and apply a correction to convert these JD values to the appropriate HJD, taking 



1 We use images processed using version S13.0 of the standard Spitzer pipeline, to avoid the additional 
noise introduced by a darkdrift correction that was applied to all subarray images beginning with version 
S14.0 of the pipeline, released in May 2006. This darkdrift correction is poorly constrained for subarray 
images dominated by a single bright star, and as a result introduces noise at a level higher than the effect 
it is meant to correct. See IRAC pipeline history available online at http://ssc.spitzer.caltech.edu/irac/ for 
more information. 
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into account Spitzer's orbital position at each point during the observations. As a check we 
repeat our analysis using apertures ranging from 3.5 — 7 pixels, and obtain consistent results 
in all cases. 

Fluxes measured at these two wavelengths show a strong correlation with the chang- 
ing position of the star on the array, at a level comparable to the dept h of the secondary 



eclipse. This effect is due to a well-documented intra- pixel sensitivity (IReach et al.l 12005 



Charbonneau et al.l 120051 ; iMorales-Calderon et al.l 120061 ) , and can be removed by fitting the 
data with a quadratic function of x and y position, where position is measured as the distance 
between the peak of the star's point spread function and the center of the pixel containing 
this peak: 



f 1 = f * ( Cl + c 2 * ( x - 14.5) + c 3 * (x - 14.5) 2 + c 4 * (y - 14.5) + c 5 * (y - 14.5) 2 ) 



(1) 



where / is the original flux from the star, f 1 is the measured flux, x and y denote the 
location of the center of the star on the array, and c\ — C5 are the five free parameters in the 
fit. We find that adding higher-order terms to this equation does not improve the fit, nor 
does adding a linear or quadratic function of time. We fit this function to the out-of-transit 
data alone and also simultaneously with the transit curve, and obtain consistent results in 
both cases. We chose to use the simultaneous fit, as it allows us to accurately estimate the 
additional uncertainty in the depth of the eclipse introduced by this correction. 

We fit the correction for the intra-pixel sensitivity of the array a nd the transit curve si- 
mult aneously to the data using a Markov Chain Monte Carlo method (lFordll2005l ; IWinn et al. 
20071 ) with 10 6 steps. We set the uncertainty on individual points equal to the standard de- 
viation of the out-of-transit data after correction for the intra-pixel variations, and remove 
outliers of 5a or more as calculated using the residuals from the best-fit light curve. We 
allow both the depth and timing of the secondary eclipse to vary independently for the 
eclipses at each of the two observed wavelengths, and take t he other parameters fo r the 
system (planetary and stellar radii, orbital period, etc.) from iKnutson et all (j2007al ). We 



calculate our transit curve using the equations from Mandel fc Agol ( 2002h for the case with 



no limb-darkening. After running the chain, we search for the point in the chain where the 
X 2 value first falls below the median of all the \ 2 values in the chain (i.e. where the code 
had first found the best-fit solution), and discard all the steps up to that point. We take the 
median of the remaining distribution as our best-fit parameter, with errors calculated as the 
symmetric range about the median containing 68% of the points in the distribution. The 
distribution of values was very close to symmetric in all cases, and there did not appear to 
be any strong correlations between variables. Figure [TJ shows the binned data with the best 
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fit to the detector effects overplotted, and Figure [2] shows the binned data once these trends 
are removed, with best-fit eclipse curves overplotted. Best-fit eclipse depths and times are 
given in Table [TJ 



2.2. 5.8 and 8.0 fim Observations (Si:As Detector) 

At longer wavelengths the flux from the star is smaller and the zodiacal background 
is larg result we chose to use a smaller 3.5 pixel aperture at these two wavelengths 

in order to minimize the noise contribution from this increased background. As a test 
we also tried using a psf fit to derive the time series in the 8 //m channel, which has the 
highest background, using the in-flight point response functions generated from calibration 
test dataB There was no improvement in the resulting time series, indicating that aperture 
photometry is still appropriate here. As a check we repeat our analysis using apertures 
ranging from 3.5 — 5 pixels and obtain a consistent signal in all cases, but with a scatter that 
increases with the radius of the photometric aperture. As before, we calculate the position 
of the star individually in each image as the position-weighted sum of the fluxes in a 7 x 7 
pixel box, and estimate the background using a Gaussian fit to a histogram of the pixels in 
the corners of the array. Fluxes in the first 10 images in each set of 64 are consistently below 
the median value for the set by as much as 10%, with the lowest values at the beginning of 
each set, so we chose to exclude the first 10 images from each set of 64 in our analysis. 

There is no intra-pixe l sensitivity at these wavelengths, but there is another well- 



documented detector effect (IKnutson et al.ll2007bl ) which causes the effective gain (and thus 



the measured flux) in individual pixels to increase over time. This effect has been referred to 
as the "detector ramp" , and has also been obser ved to occur in the IRS an d MIPS 24 /jm ar- 
rays, which are made from the same material JPeming et al.lbood . liooS l. The size of this 



effect depends on the illumination level of the individual pixel; pixels with high (>250 MJy 
Sr~l in the 8 /xm channel) will converge to a constant value within the first hour of obser- 
vations, while lower-illumination pixels will show a linear increase in the measured flux over 
time, with a slope that varies inversely with the logarithm of the illumination level. 

This effect is important for two reasons. First, it means that the observed 3% linear 
increase in the measured background flux at 8 /im over the period of the observations is most 
likely not the result of a real change in the zodiacal background, but is instead another exam- 
ple of this detector ramp. Although the increased noise and smaller size of the background 
at 5.8 /im obscures this effect, there appears to be a similar upward trend. Thus, rather 



2 Available at |http://ssc.spitzer.caltech.edu/irac7 psf.html 
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than calculating the background in each image individually and subtracting that value, we 
subtract a constant background of 3.21 MJy Sr~l per pixel from all the 8 /im images, and 
0.41 MJy Sr _1 per pixel from all 5.8 /im images. This background is calculated as the median 
background value during the last 2.5 hours of observation, when presumably the background 
is closest to its true value. We note that this choice has a negligible effect on our final eclipse 
depths, as the background constitutes only 0.2% and 2.3% of the signal in our aperture at 
5.8 and 8.0 /im, respectively. If we subtract no background at all, the eclipse depths we 
obtain are still well within the la uncertainties of our final quoted values. 

This effect also produces a 0.5% increase in the measured flux from the star at 8.0 /im over 
the period of these observations (see Figured]), and a much smaller (< 0.1%) increase at 
5.6 /im. Unlike the detector ramp at low illumination levels, the ramp for higher illumina- 
tions has an asymptotic shape, with a steeper rise in the first 30 minutes of observations. We 
discard the first 30 minutes of data in both the 5.8 and 8.0 /tm channels and fit the remaining 
binned time series from our 3.5 pixel aperture with a qua dratic function o f ln(dt) , where dt 



is the change in time from the start of the observations. iKnutson et al.l (l2007bl ) used this 
same functional form (with additional degrees of freedom) to describe the detector ramp in 
the 8 /im channel over a period of 33 hours of continuous observations of HD 189733, and 
it accurately c a ptures the behavior of this ramp for a range of illumination levels. Unlike 
Knutson et al. (b007bh we do not attempt to correct each of the pixels in the images indi- 



vidually for this ramp; this is not neccessary for our analysis, and the lower fluxes, shorter 
time frame (8 hours instead of 33) and reduced cadence of our data (from cycling between 
the four detectors) make it difficult to characterize this effect accurately at the pixel level. 
Instead we assume the detector ramp in the binned flux in our 3.5 pixel aperture will have a 
shape similar to the ramp for individual pixels and fit the binned time series in each channel 
with a quadratic function of ln(dt), which produces a good fit to the observed trends. 

We fit both the quadratic function of ln(dt) and the transit curve to the data simul- 
taneously using a Markov Chain Monte Carlo method as described in §2.11 As before, the 
distribution of values was very close to symmetric in all cases, and there did not appear to 
be any strong correlations between variables. Best-fit eclipse depths and times from these 
fits are given in Table [TJ and the time series both before and after correcting for detector 
effects are shown in Figures [T] and [21 respectively. As a check we repeated these fits using 
both a linear function in time and a linear function of ln(dt), and found that the eclipse 
depth in both cases varied by less than la from our quoted value. 
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Discussion 



We ultimately achieve noise levels in each of the four bandpasses that are 1.5, 1.9, 2.0, 
and 1.6 times higher than the predicted photon noise at 3.6, 4.5, 5.8, and 8.0 /xm, respectively. 
The additional noise is most likely the result of the jitter introduced by the need to repoint 
the telescope every 3.5 minutes in order to switch between bandpasses. Although the eclipse 
is detected to a high degree of significance in all four bandpasses, it is worth noting that 
we ultimately achieve a precision at 4.5 and 8.0 tim comparab le to that of the measured 
secondary eclipse depths for TrES-1 (ICharbonneau et al.l 120051 1. even though this star is 
significantly fainter (K = 9.82 vs K = 6.32 for HD 209458). This is partially explained by 
the reduced cadence of the HD 209458 observations, which had 15% of the total effective 
integration time per band for the TrES-1 observations, but the f requent repointings appear 
to have contributed additional noise as well. We also note that ICharbonneau et al.l (120051 ) 
did not include the uncertainties contributed by their fits to the trends in the out-of-transit 
data in their error estimates. 

If we combine our estimates of the eclipse timing in each of the four bandpasses, we find 
that the center of the eclipse occ urs at 2453702 5251 ± 0.0012 HJD. This is 4.9 ±1.7 minutes 
or 2.9a earlier than predicted (IKnutson et al.l l2007al ). and this is without accounting for 
the addition al 50 s delay in the predicted time due to the light travel time in the system 
( lLoebll2005l ). In this case the uncertainty in the predicted time is negligible compared to the 
uncertainty in our measurement. 

It is possible that we have under-estimated the uncertainties in our estimates for the 
timing of the eclipse. Because the ingress and egress occur on relatively short time scales, 
fitted values for the eclipse times are particularly sensitive to the presence of correlated noise 
in the time series, which is not included in the error bars from the Markov fits described 
above. To estimate the effect of correlated nois e on our best - fit ecl ipse depths and times, we 
use the "prayer-bead" method as described in iGillon et al.l (120071 ). In this method we take 
the time series of the residuals from our best-fit solution for each eclipse, shift the residuals 
forward in time with the points at the end of the time series wrapping back around to 
the beginning, add the best-fit solution back in, and fit the new timeseries for the full set of 
parameters. This process is iterated until the residuals have been shifted back around to their 
original positions in the timeseries. The variance in the resulting set of values for the best-fit 
eclipse depths and transit times gives the error values for each of the parameters. In this 
case we set our uncertainties equal to the range in values containing 68% of the points in the 
distribution in a symmetric range about the median for a given parameter, as we did before 
with our Markov fits in §2.11 The resulting uncertainties are comparable to the uncertainties 
we obtained from our Markov fits, with some larger and some smaller values. For the eclipse 
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depths, we find uncertainties of [0.0009, 0.00025, 0.00072, 0.00020] at [3.6, 4.5, 5.8, 8.0] fj,m. 
For the best-fit eclipse times, we find uncertainties of [4.0, 2.9, 5.4, 1.9] minutes. From 
this we conclude that correlated noise is not a significant source of uncertainty in the data. 
Significantly, we note that the uncertainty in the best-fit time for the 4.5 /im eclipse, which 
occurs 12.6 minutes earlier than predicted, is smaller than the original uncertainty from our 
Markov fit. We elect to use the uncertainties from our Markov fits as the final uncertainties 
for our parameters, as we feel that this is the better method in the case where correlated 
noise in the data is minimal, as it appears to be. In either case the differences between the 
two methods for calculating uncertainties are minor and do not affect our conclusions. 

Although the combined best-fit time for the four eclipses appears to occur 4.9 minutes 
or 2.9cr earlier than predicted, we do not believe that this is a convincing detection of a non- 
zero orbital eccentricity. Because we observe the eclipse simultaneously at four wavelengths, 
we would expect to see similar timing offsets in all four channels if the shift was the result 
of a non-zero eccentricity. Instead, there appears to be a larger, marginally significant 
(3.6a) offset in the 4.5 fim channel, while the other three channels are effectively consistent 
with zero offset (see Table [Q . A more plausible explanation would be an apparent timing 
shift caused by a color-dependent non-uniform brightness distribution on the surface of 
the planet, which would alter the shape of t he ingress and egress relative to the shap e 
expected for a uniform brightness distribution (jWilliams et al.l 12006c iRauscher et al.l 120071 ). 
It is reasonable to expect that this brightness distribution might vary with wavelength, 
causing different apparent timing shifts in each of the four channels, as different wavelengths 
probe different depths in the atmosphere. The signal-to-noise of our data is not high enough 
to distinguish the changes in the shape of the ingress and egress that would signal such a 
non-uniform brightness distribution, but continuous, higher-cade nce observations similar to 
the 8 jum observations of HD 189733b by lKnutson et al.l (j2007bl ) might provide a definitive 
answer to this question in the future. 

When we compare the measured eclipse depths at 3.6, 4.5, 5.8, and 8.0 /im to the 
predictions of a theoretical model for this planet (see Figure [3]), it is immediately clear that 
the shape of the observed spectrum differs significantly from the predicted values. The data 
in Fig. [3] show a peak in flux centered around the 5.8 /im bandpass and indicate that the 
flux in the 4.5 fim bandpass exceeds that in the 3.6 /im bandpass. As Fig. [3] suggests, 
the effective photospheres of the 4.5 and 5.8 /im features are at rather high temperatures 
in the atmosphere and the correspo nding effective temperature of the 3.6 /im flux is at a 



lower temperature. Prey ious theory (IBurrows et al.ll2005c iFortney et al.ll2005t iBarman et al. 



2005c ISeager et al.ll2005l ) had suggested that there would be a trough between the 3.6 and 
8.0 fxm IRAC bandpasses due to a water absorption feature and that the flux at 3.6 //m would 
either exceed that at 4.5 fim or would be comparable to it. However, the new IRAC data 
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force us to conclude that there is a thermal inversion in the atmosphere of HD 209458b at 
altitude and that we are indeed seeing water in the 4 — 8 //m gap, but in emission. In fact, 
we now expect the water absorption features of the older default theory to be flipped into 
emission features throughout the entire n ear- to mid-infrared wav elength range. Therefore, 
we find that, contrary to the conclusion of iRichardson et al.l (120071 ). the flatness or slight rise 
of their IRS spectrum near ~7.8 /im in fact supports the presence of abundant atmospheric 
water, because this spectral region is at the edge of a strong water band in emission. A 
temperature inversion might a l so na tu rally explain the two emission features tentatively 
identified by IRichardson et al.l (120071 ). iBurrows et al.l (j2007bl ) explore the theoretical and 
model consequences of these new IRAC data in more depth. 

The idea that the spectrum of a strong ly irradiated extrasol ar gi ant planet could man i- 
fest water emission features was presaged by lHubeny et al.l (120031 ) and lBurrows et al.l (120061 ). 
who discussed a bifurcation in the atmosphere solution, b ut we still do not know t he nature 
of the stratospheric absorber responsible for this heating. iRowe et al.l (120061 . 120071 ) report a 
value of 0.04 ± 0.04 for the geometric albedo of the planet at visible wavelengths, indicating 
that the stra tospheric abso r ber m ust absorb much more in the optical than the infrared, as 
predicted by iHubeny et al.l (120031). This eliminates m any types of clouds, which tend to be 
reflective at optical wavelengths. ICowan et al.l (120071 ) placed a 2a upper limit of 0.0015 on 
the size of the phase variation for this planet at 8 /tm; when compared to the eclipse depth 
of 0.00240 ± 0.00026 described in this paper, this indicates that the flux from the night side 
must be at least 60% of the flux from the day side, which would indicate relatively efficient 
circulation in the atmosphere at the level of the 8 /mi photosphere. Because high-altitude 
clouds or other opaque layers would shift the location of the 8 jum photosphere, this provides 
additional constraints on the nature of the unknown stratospheric absorber. 



4. Conclusions 

We estimate the secondary eclipse depth for the transiting planet HD 209458b at 3.6, 4.5, 
5.8, and 8.0 /im. These observations provide a useful complement to previous observations 
of this system, which were limited to wavelengths longer than 8.0 /mi. In contrast to the 
results for longer wavelengths, we find that the planet's emission at shorter wavelengths 
is clearly inconsistent with the predictions for a standard cloudless atmosphere model. We 
suggest an alternative explanation, in which a temperature inversion in the upper atmosphere 
produces water emission features at 4.5 and 5.8 /im. Although the cause of this inversion 
layer is unknown, it is suggestive that this planet also falls into a class of hot Jupiters which 
appear to have radii that are significantly larger than predicted by standard models for an 



- 11 - 



irradiated gas giant. These two characteristics may or may not be related; the obvious next 
test would be to extend these observations to a much larger sample of planets to determine 
which, if any, of the other 21 known transiting planets show similar tempe rature inversions 



The v ery bright 8.0 fim flux from the core-dominated planet HD 149026b ([Harrington et al. 



2007 ) indicates that this planet may also have a temperature inversion, as predicted by 



Fortney et al.l (j2006bl ); this result should be confirmed by Spitzer observations at additional 



wavelengths in the near future. 

Given the imminent depletion of Spitzer's cryogen, at which point only the 3.6 and 
4.5 /im channels will be functioning, it is worth noting that the best evidence for the tem- 
perature inversion in HD 209458b's atmosphere comes from observations in these two shorter- 
wavelength channels. Spitzer provides the optimal platform for this type of measurement, 
and it would be relatively straightforward to survey all of the known bright transiting sys- 
tems as part of the post-cryogenic mission^]. Such a survey has the potential to provide 
a definitive answer to the question of what properties of the planet or parent star lead to 
these temperature inversions, and perhaps shed light on the nature of the clouds or other 
upper-atmosphere optical absorbers that are needed to produce the temperature inversion. 
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Table 1. Best-Fit Eclipse Depths and Times 



A (/im) 


Relative Eclipse Depth 


Center of Transit (HJD) 


O-C (min.) b 


3.6 


0.00094 ± 0.00009 


2453702.5244 ± 0.0024 


-5.9 ±3.4 


4.5 


0.00213 ± 0.00015 


2453702.5198 ± 0.0024 


-12.6 ±3.5 


5.8 


0.00301 ± 0.00043 


2453702.5251 ± 0.0026 


-4.9 ±3.8 


8.0 


0.00240 ± 0.00026 


2453702.5299 ± 0.0022 


2.0 ±3.1 


24 a 


0.0026 ± 0.00046 


2453346.5278 ± 0.0049 


-1.6 ± 7.1 



HPemine et"al] j2005h 

b Predicted transit time from lKnutson et al. | l l2007al'l 
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Fig. 1. — Secondary eclipse of HD 209458b on UT 2005 Nov. 27, observed at (from top to 
bottom) 3.6, 4.5, 5.8, and 8.0 fim, binned in 7 minute intervals and normalized to one. The 
overplotted curves show the best-fit corrections for detector effects (see §2.11 and §2.21) . 
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Fig. 2. — Secondary eclipse of HD 209458b on UT 2005 Nov. 27, observed at (from top 
to bottom) 3.6, 4.5, 5.8, and 8.0 fim, with best-fit eclipse curves overplotted. Data has 
been normalized to remove detector effects (see discussion in §2.11 and §2.2|) . and binned in 
7 minute intervals. 
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Fig. 3. — Predicted emission spectrum for HD 209458b from lBurrows et al.l (120061 ). Squares 
show this spectrum integrated over the Spitzer bandpasses (response functions for these 
bandpasses are shown at the bottom of plot (dotted lines), scaled by a factor of 0.001), 
and circles are the estima t ed ec lipse depths in these bandpasses. The 24 fim point is 
taken from iDeming et al.l (120051 ). The two dashed lines give the planet-star flux ra- 
tio for the case where the planet is a perfect blackbody with a temperature of either 
1500 K or 1900 K and the spectrum of the star is calculated from a model (available 
at http://kurucz.harvard.edu/stars/hd209458). Note that w e have chosen not to plot the 



7.5 — 12.2 fim spectrum measured by iRichardson et al.l (120071 ): this is because the spectrum 
was scaled to match a preliminary value for our 8.0 ^m eclipse depth, and thus does not 
contain independent information about the absolute strength of the emission from the planet 
at these wavelengths. 



